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The ferroelectric trends of rutile (Ti02) in bulk and at the (110) surface are investigated by 
means of ah initio density functional theory. We discuss the underlying mechanism of the incipient 
ferroelectric behavior of rutile in terms of Born effective charges, which we decompose in individual 
contributions by means of maximally localized Wannier functions. We show that a ferroelectric 
phase can be stabilized for a variety of different lattice distortions which all enlarge the shortest 
Ti-0 distance, even if the longer apical Ti-0 bond is simultaneously shortened. At the (110) surface, 
the ferroelectric trends are modified compared to the bulk, but nevertheless ferroelectric phases with 
large polarization even in the topmost surface layer can be stabilized by uniaxial strain. 



I. INTRODUCTION 

Ti02 in the rutile structure has been classified as an 
incipient ferroelectric material, based on the experimen- 
tal observation that the dielectric permittivity is excep- 
tionally large and increases with decreasing tempera- 
ture, although no ferroelectric transition occurs even for 
T = K.^ ^ The realization of ferroelectric Ti02, or even 
the possibility to efficiently tune its dielectric constant, 
would be of great technological importance, since Ti02 
is a cheap material with well optimized processing meth- 
ods, which is used for a variety of applications such as 
optical coating,-"^ solar cells, ^ or sensors.-^ Many of these 
applications are based on properties involving the rutile 
(llO)-surface, which is the technologically most relevant 
and therefore best-studied rutile surface .I^^KI^ 

Previous theoretical ah initio studies have shown that 
in the corresponding bulk material the anomalous be- 
havior of the dielectric constant is a consequence of a 
low-frequency polar (A2n) phonon mode,!^ which softens 
under negative pressure or uniaxial tensile strain along 
the (00 1) -axis .l^^^^ Liu et al. classified the phase transi- 
tion between the paraelectric and ferroelectric phases to 
be of second order. Since the resulting magnitude of the 
ionic off-centering as well as the electric polarization are 
in the same range as for prototypical ferroelectrics such 
as PbTiOa or BaTiOs, an enforced ferroelectric transi- 
tion would lead to a technologically relevant ferroelectric 
material. Furthermore, Mitev et al. found a large re- 
gion within the Brillouin zone with a low lying phonon 
branch which may be linked to the structural relaxation 
of the rutile (110) surface. These authors also pointed 
out a softening of this mode under uniaxial strain in both 
[001] and [110] directions, and found that the ferroelec- 
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trie A2n-niode is quite insensitive towards [110] -oriented 
strain. In contrast, we have shown more recently that 
the A2n-mode indeed also softens under [110] strain, but 
that a different polar mode, an mode polarized along 
[110], becomes even more unstable than the A2n mode 
under this strain orientation.'^ Under uniaxial tensile 
[110] strain, we have thus found a second order phase 
transition to a ferroelectric phase with polarization along 
[110].^^ In agreement with these results, a softening both 
of the A2n mode and of the lowest E^^ mode has also been 
found in Ref.|15|for tensile uniaxial [001] strain and for an 
isotropic lattice expansion within the (100)/(010) plane. 
Since the strain dependence of the two polar modes is dif- 
ferent, this opens up the interesting possibility to "strain 
engineer" the polarization direction. 

It is well known that ferroelectric or dielectric prop- 
erties can be significantly modified at surfaces in com- 
parison with the corresponding bulk material. First of 
all, a polarization along the surface normal induces a 
depolarizing field, and is therefore often suppressed in 
thin films (see, e.g., Ref. [19]). Second, the reduced co- 
ordination of the surface atoms may prevent or enhance 
ferroelectric trends as the dipolar interaction and 
the short-range repulsion are both modified. Finally, 
the surface-induced atomic relaxations can also modify 
the bonding and thus the ferroelectric properties. For 
example, it has been discussed whether the surface re- 
laxation and the unsaturated bonds at the surface are 
able to destabilize the paraelectric state, leading to a 
finite electric polarizatio n in th e related incipient ferro- 
electric material SrTiOs.l ^ ' ^^ ' ^^ ' We also note that a fer- 
roelectric transition has been observed experimentally in 
strained SrTiOa films. Therefore, an obvious approach 
to achieve a ferroelectric transition in rutile is the growth 
of strained thin Ti02 films. 

In view of this, it is essential to test the modifications 
of the ferroelectric properties of the strained Ti02 (110) 
surface, as considerable differences may appear compared 
to the bulk case. Nevertheless, to the best of our knowl- 



2 



edge no such investigation of the ferroelectric trends of 
this Ti02 surface has been undertaken until now. 

In the present paper we first present a systematic com- 
parison of the effects of different lattice modifications on 
the ferroelectric properties of bulk rutile Ti02. We show 
that each lattice modification which affects the nearest 
neighbor Ti-0 distances within the crystal, can stabilize 
a ferroelectric phase polarized along (001). Additionally, 
we discuss the fundamental mechanisms for the incipient 
ferroelectric behavior of rutile Ti02 in the framework 
of Born effective charges. We then present the ab ini- 
tio investigation of the ferroelectric trends of the Ti02 
rutile surface. Although we find strong surface-induced 
relaxations and modifications of the electronic structure, 
which have a considerable influence on the ferroelectric 
properties, the strong ferroelectric tendencies of the bulk 
material are essentially retained at the surface. In addi- 
tion, two different ferroelectric states with polarization 
within the surface plane can be stabilized for strained 
Ti02 films. Here, local dipoles both along (001) and 
along (110) can be stabilized within the surface at 2 % 
uniaxial strain. Furthermore the magnitude of the polar- 
ization increases as the strain is enlarged to 5 %. 

This paper is organized as follows. In Sec. [IT] we 
present our computational methods and all related tech- 
nical details. In Sec. lIIII we discuss the ferroelectric trends 
of bulk rutile in some detail. This is essential in order 
to make reliable predictions on the ferroelectric trends of 
the surface discussed in the later sections. For this pur- 
pose, the ferroelectric trends under different lattice mod- 
ifications and a discussion of the underlying mechanisms 
will be presented in Sec. |III.2] and Sec. |III.3l respectively. 
We outline the ferroelectric trends of the rutile (110) sur- 
face in Sec. IIVI both in free films and for films which are 



clamped to an idealized substrate (Sec. IV.3). Finally, we 
show that different ferroelectric distortions can be sta- 
bilized in the surface under tensile uniaxial strain and 



discuss the resulting displacement patterns in Sec. IV. 4 



Conclusions and outlook on future work are summarized 
in Sec. m 



II. COMPUTATIONAL DETAILS 



dielectric properties compared to an enlarged valence ba- 
sis set and thus stick to the minimal valence electron ba- 
sis. A planewave cutoff of 500 eV, an energy convergence 
better than 10"'^ eV, and at least 6 x 6 x 6 (5 x 13 x 3) 
k-points constructed with the Monkhorst-Pack^^ scheme 
for bulk (films) guarantee sufficient accuracy of our calcu- 
lations. A Gaussian smearing of 0.1 eV is used for most 
calculations in order to improve the convergence with 
respect to the number of k-points. This smearing is re- 
duced to 0.05 eV for the determination of the bulk Born 
effective charges, and the linear tetrahedron method is 
used for the determination of the harmonic frequencies. 
Atomic positions are optimized until the residual forces 
are smaller than 0.001 eV/A (0.005 eV/A) for bulk (sur- 
face) calculations. For the surface calculations at least 
10 A of vacuum is used between the periodic images of 
the films. 

The Plane- Wave Self-Consistent Field (PWscf) code, 
which is part of the QuantumEspresso package, in com- 
bination with Wannier90,^ is used to construct maxi- 
mally localized Wannier functions (MLWFs)^^ for the or- 
bital decomposition of the Born effective charges of bulk 
rutile Ti02. For this purpose, a F-centered equidistant 
14 X 14 X 20 /c-point grid, a plane-wave cutoff of 35 Ry 
(^ 476 eV) for the wave- functions (420 Ry for the charge 
density), and an energy convergence of 10~^ Ry ensure 
sufficient accuracy. The Ti 35 and 3p semi-core states are 
treated as valence electrons in the PWscf calculations. 
For the MLWF decomposition of the surface, the /c-point 
grid is reduced to 6 x 12 x 2. MLWFs are constructed sep- 
arately for valence and semi-core states, based on appro- 
priate initial projections, and the total spread is always 
convere; ed better than 10"^° (10"^ A^) per Wannier 
function for the bulk (surface). 



III. BULK RUTILE 

In order to provide a well-defined reference for the in- 
terpretation of our results for the Ti02 rutile (110) sur- 
face in Sec.[TVl we first discuss results for the correspond- 
ing bulk system. 



The electronic structure of Ti02 is calculated using 
first-principles density functional theory. Most calcula- 
tions are performed by employing the Vienna Ab Ini- 
tio Simulation Package (VASP 5.2.2)^^ based on pro- 
jector augmented wave pseudo-potentials.^^ Born effec- 
tive charges and the dielectric permittivity are calcu- 
lated within VASP using density functional perturba- 
tion theory. The generalized gradient approximation 
(GGA) in the formulation of Perdew, Burke and Ernzer- 
hof (PBE)^^ is used for the exchange correlation energy, 
with 4:s3d {2s2p) electrons treated as valence for Ti (O). 
Although small inaccuracies may result from the use of 
such a soft Ti pseudopotential,^ we could not find qual- 
itative differences in the description of the electronic and 



III.l. Structural properties 

The atomic structure of Ti02 in its rutile morphol- 
ogy is tetragonal with PA2mnm symmetry, see Fig. [l] 
Besides the lattice constant a and the tetragonal c/a- 
ratio, an internal parameter u defines the Ti-0 dis- 
tances. Each Ti-atom is 6-fold coordinated by 0-atoms 
with 4 short equatorial Ti-0 bonds (Ti-Oeg), deq = 
|(22x — 1, 2ii — 1, c/a)| • a/2, and two long apical Ti-0 
bonds (Ti-Oap), dap = \{u,u,0) \ - a. 

It has been shown that the ferroelectric trends 
of bulk rutile de pend critically on the lattice 
parameters.l^^^I^HIS^ Because of this, a correct de- 
scription of the atomic structure is essential. Our 
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FIG. 1. (Color online) Atomic structure of rutile. Ti (blue, 
light), O (red, dark), (a) Primitive unit cell. Arrows mark the 
atomic displacement of the phonon mode. The atomic 
structure can be decomposed into O-Tis units, see red trian- 
gle, (b) ^ X 1 X ^ cell in [110], [001], [110] coordinates. 
Ti(l)/Ti(2) /0(l)/0(2) mark symmetry-inequivalent atomic 
classes under [110] strain. Arrows mark the polar atomic dis- 
placements under [110] strain, see Ref. 1181 



TABLE II. Modification of the lattice parameters under ex- 
ternal pressure P. Our values obtained with GGA-potentials 
are opposed to values obtained with LDA potentials from 
Ref. [141 £^nd experimental data is given for comparison. For 
each quantity the calculated pressure dependence has been 
fitted by /(P) = bo + b ■ P . The values for zero pressure, 6o 
(left), and the slopes 6, in units 10^ GPa~^ for u and c/a, in 
units of 10^ GPa~^A otherwise (right) are listed. 







bo 






b 




ldaE^ 


Here 


Exp.Eil 


LDAE^ Here Exp.E^ 


a [A] 


4.55 


4.66 


4.59 


-7.73 


-9.1 -8.62 


c[A] 


2.92 


2.97 


2.96 


-2.41 


-2.2 -2.14 


U 


0.304 


0.305 




-0.14 


-0.2 


c/a 


0.64 


0.64 


0.64 


0.56 


0.8 0.79 


TiOeg[A] 


1.93 


1.97 




-1.70 


-1.5 


TiOap[A] 


1.96 


2.01 




-4.27 


-5.4 



TA BLE I . Bulk lattice parameters and Born charges Z* (see 
Sec. HI. 3 for further details) of rutile Ti02. Results obtained 
in the present work (top rows) are compared to literature data 
(bottom rows). 



a [A] c/a u 



(001) 



GG^ 


a 


4.664 0.637 0.305 7.60 5.36 


8.15 


GGA 


b 


7.25 - 


7.73 



LDA^ 4.572 0.644 0.304 7.49 5.43 7.77 

LDA^^ 4.535 0.641 0.303 7.27 5.25 7.32 

GGA^^ 4.637 0.638 0.305 7.22 5.18 7.52 

exp[ ^ 4.587 0.644 0.305 - 

^ VASP 

^ PWscf; structural parameters fitted to VASP results. 
^ 4.2 K, Ref.l2] 



calculated lattice parameters, see Table [l[ and the 
calculated pressure dependence of these quantities, 
see Tab. [llj fit better to experimental results than 
calculations based on the commonly used local density 
approximation (LDA). Nevertheless, the LDA error 
is over corrected, leading to an overestimation of the 
lattice constant by 1.6% compared to low tempera- 
ture experimental results, which may also lead to a 
quantitative overestimation of the ferroelectric trends. 
We note that, while an artificial destabilization of 
the paraelectric phase for GGA potentials has been 
found in Ref. 13^ we could not reproduce this in our 
calculations, in agreement with Ref. |35| see Tab. Ill 
These diflFerences can most likely be attributed to the 
diflFerent pseudopotentials and technical details such as 
the used basis functions, see also the related discussion 
for the incipient ferroelectric material SrTiOs.'^ 




FIG. 2. (Color online) Stiffness s of the mode for differ- 
ent modifications of the equilibrium crystal structure: Exter- 
nal pressure (blue, stars); i^- variation (red, diamonds); [001] 
strain (cyan, squares); [110] strain (green, triangles); isotropic 
lattice expansion (black, crosses). Left inset: Energy differ- 
ence (meV/2 f.u.) for a static atomic displacement at the 
equilibrium lattice constants (black, diamonds) and under 1 % 
tensile [001] strain (red, squares). /\z is the relative Ti-0 dis- 
placement in [001] direction. Right inset: Ratio of the Ti-Oap 
and Ti-Oeg distances under the imposed lattice modifications. 



III. 2. Polar phonon modes 

In the following, the relative stability between the 
paraelectric state and a ferroelectric state with polariza- 
tion along [001] will be discussed in the framework of the 
frozen phonon approach, cf. Refs. 13] [HI [161 and (361 For 
this purpose the atoms are successively displaced along 
the eigenvector Y of the A2u phonon mode and the total 
energy is monitored. 

The eigenvector of the A2n mode is fully determined by 
symmetry as it is the only polar mode along [001],^ and 
a gradual shift of amplitude A along the eigenvector Y 
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TABLE III. Coefficients of total energy fits, Eq.[T] for atomic 
displacements along the A2u mode and harmonic frequencies 
cjo in comparison with literature data. Top part: external 
pressure; Center: tensile [001] strain; Bottom part: tensile 
[110] strain. 





Equation 1 


Harmonic frequencies 




s 




PBE LDA^i^ PBeP LDA'^ 




(eV/A^) 


(eV/A^) 




cJo(cm~^) 


Eq. 


0.51 


21.89 


91 


150 i86 170 


-2 GPa 


0.14 


21.36 


47 


100 - 150 


-5 GPa 


-0.58 


21.71 


i 97 




1 % 


-0.50 


20.74 


i90 


130 


2 % 


-1.36 


17.85 


il48 


100 


5 % 


-4.00 


15.29 


i254 


il25 


2 % 


-0.76 




illO 




5 % 


-2.30 




il92 





is thus given by the displacements A - (0, 0, 1) for Ti and 
A - (0, 0, —rriTi/^mo) for 0-atoms, respectively. Here, 
represent the Ti and O masses. 

While the paraelectric state is stable at the equilib- 
rium lattice constants, the energy surface flattens around 
^4 = for tensile [001] strain, or under an isotropic vol- 
ume expansion. In both cases the A2u mode softens pro- 
portional to the imposed lattice modification, and two 
ferroelectric minima appear at a critical value of the lat- 
tice modification,^^ see left inset in Fig.|2j Although 
so far no such transition has been observed experimen- 
tally, a considerable modification of the A2u mode fre- 
quency under external pressure has been measured.^ 

In order to discuss these ferroelectric trends in more 
detail, we expand the energy change per primitive cell 
during the polar displacement, A£^(A), up to 4^^ order 
with respect to lA • y|. 



Y 



kA^ 



Y 



(1) 



For a qualitative discussion we use the stiffness s, which 
determines the curvature of the energy at the paraelectric 
state. 

The stiffness is positive for a stable paraelectric state 
and changes its sign as the ferroelectric state becomes 
the energetic ground state. Figure [2] depicts the vari- 
ation of the stiffness as a function of the Ti-Ogg bond 
length for the following modifications of the equilibrium 
crystal structure: {%) external pressure, {%%) isotropic lat- 
tice expansion, {iii) [001] strain, (iv) [110] strain, and 
(v) variation of the ratio between TiOgg and TiO^p bond 
lengths at constant volume by a variation of the internal 
parameter u. Here, for {i)^ all lattice parameters have 
been optimized at a given pressure, whereas neither the 
lattice vectors nor the internal parameter u have been 
optimized for all other cases. 

It can be noticed that we obtain a stable ferroelec- 
tric state polarized along [001] for each tested lattice 



modification which enlarges the Ti-Oeg bond, see Fig. [2j 
Even an artificial modification of the internal parame- 
ter li, which enlarges the Ti-Ogg bonds and reduces the 
Ti-Oap distances without lattice expansion, stabilizes the 
ferroelectric phase. 

For the case of [001] strain, it has been discussed in 
Ref. 114, that the large short-range repulsion, which ap- 
pears under a relative displacement of the two sublat- 
tices, prevents a ferroelectric transition for the equilib- 
rium structure, but that this contribution to the total 
energy decreases with increasing strain. Our system- 
atic comparison of different lattice modifications further 
shows that it is sufficient to increase the Ti-Ogq distance 
in order to stabilize the ferroelectric phase. Most likely, 
the relative displacements of these ions contribute most 
to the short-range repulsion as they correspond to the 
shortest Ti-0 distance. 

Even though all the various lattice modifications desta- 
bilize the paraelectric phase, the slopes of 5 as a function 
of the Ti-Ogg bond length differ considerably for the dif- 
ferent cases. Here, the differences in the Ti-Oap distances 
result in different contributions of these atomic pairs to 
the short-range repulsion. Additionally, the dipolar in- 
teraction is reduced as the ratio of the nearest and next- 
nearest Ti-0 distances changes, see Sec. |IIL3[ As the 
dipole-dipole interaction is the second large contribu- 
tion to the to tal en ergy, that is affected by ferroelectric 
displacement s,'^^^ its reduction reduces the ferroelectric 
trends considerably. Because of this, the slopes of s are 
correlated with the ratio of the two classes of Ti-0 bond 
lengths, see right inset of Fig. [2] On the one hand, the 
configuration for which the atomic structure has been op- 
timized under external pressure, shows the steepest s in 
Fig. |2j as an increase of the ratio of the TiO^p and TiOeg 
distances is superimposed to the lattice expansion, see 
Tab. [n[ On the other hand, s is notably flatter if the ra- 
tio of the two Ti-0 distances decreases, e.g., under [001] 
strain or especially for a modification of which induces 
the largest decrease of the bond length ratio. 

For a quantitative comparison of the obtained mode 
stiffness, the harmonic frequencies cjo can be obtained 
from Eq. [l] 



- 2mTi + 4 • mo 



rriTi 
2 • mo 



(2) 



In this way, we obtain a frequency of the A2u mode at the 
equilibrium structure, see Tab. |III[ which is within the 
spread of different theoretical estimates and below the 
experimental value of 173 cm~^.^ See also the discussion 
in Sec. III.1| Apart from these quantitative differences, 
similar qualitative trends for the mode softening under 
a lattice modification (pressure and [001] strain) are ob- 
tained within the different theoretical studies. 

In summary, a ferroelectric phase polarized along [001] 
can be stabilized for each tested expansion of the Ti- 
Oeq bond. We will see in Sec. |IV.4| that this result can 



qualitatively be transferred to the ferroelectric trends of 
the surface. 
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TABLE IV. Modification of tfie internal parameter u, polar 
displacements z of the Ti and O ions in A, and the energy 
difference relative to the optimized paraelectric state in 
meV/atom for different imposed tensile strains. Ti(l)/0(1) 
corresponds to Ti-Oap bonds along [110], Ti(2)/0(2) to Ti- 
Oap bonds along [110], see Fig. WCb). Upper part: Ferro- 
electric state polarized along [00 ij; Lower part: Ferroelectric 
state polarized along [110]. 



Pol. 


Strain 


^[110] '^^[110] 


ZTi(l) ZTi{2) 


ZO{l) Zo{2) 


AE 




[001]+2% 


0.305 


0.05 


-0.05 


-1 


[001] 


[001]+5% 


0.303 


0.08 


-0.13 


-11 


[iio]+2% 


0.305 


0.08 0.01 


-0.02 


-1 




[iio]+5% 


0.035 0.036 


0.14 0.02 


-0.04 -0.03 


-4 




[001]+2% 


0.306 


-0.02 -0.04 


0.03 0.02 


-0.1 


[110] 


[001]+5% 


0.307 0.306 


-0.04 -0.07 


0.07 0.04 


-2 


[110]+2% 


0.307 0.305 


-0.04 -0.10 


0.04 0.03 


-2 




[iio]+5% 


0.311 0.307 


-0.05 -0.19 


0.09 0.03 


-16 



TABLE V. Contributions of individual MLWFs (see Fig. [s] 
and total Z^^i. Oeq (Oap) correspond to MLWFs situated at 
equatorial (apical) oxygen atoms relative to the displaced Ti. 
Rows denoted "Eq." are calculated for the equilibrium struc- 
ture, rows denoted "i^+3 %" are calculated for 3 % increased 
internal parameter u. Upper part: ZJqq^j; Lower part: Z*|. In 
each case the mean values for the corresponding orbitals are 
given. 









Oeq 






Oap 










P± 


Spiq 


SPap 


P± 


Splq 


SPap 




^[001] 


Eq. 


0.48 


0.30 


0.04 


-0.01 


-0.01 


-0.10 


3.96 




1^+3% 


0.39 


0.29 


0.00 


-0.03 


-0.02 


-0.08 


3.26 




Eq. 


-0.36 


-0.04 


0.06 


0.48 


0.23 


0.29 


3.38 




1^+3% 


-0.23 


-0.08 


0.05 


0.53 


0.27 


0.31 


4.20 



III. 3. Born charges 



In the remainder of this sub-section we now discuss the 
actual displacement patterns of the atoms under differ- 
ent values of strain, which will be used as reference for 
the atomic displacements within the strained surfaces in 
Sec. |IV.4[ For this purpose, we strain the bulk structure 
along the [001], respectively [110], direction and relax all 
atomic positions, whereas we do not consider a further re- 
laxation of the lattice vectors perpendicular to the strain 
direction. If the atomic positions are optimized within 
the metastable paraelectric state, the internal parameter 
u slightly increases with strain (<1 % for imposed strain 
of up to 5 %), i.e., the ratio between the Ti-O^p and 
Ti-Oeg distances increases, as the Ti-Ogq distances are 
less sensitive towards strain or lattice expansion. This 
relaxation of the atomic positions is in agreement with 
previous work, e.g., Ref. 1141 

If the paraelectric symmetry of the structure is lifted, 
we obtain polar displacements of the atoms, in agreement 
to the discussed softening of the polar modes under ten- 
sile strain, see Table [TV| Polar displacements both along 
[001] and [110] directions lower the energy relative to the 
paraelectric state under 2 % and 5 % tensile strain, and 
the cases with polar displacements parallel to the strain 
direction are energetically most favorable for both strain 
orientations, cf. Ref 15, 

We note that the amplitudes of the polar displace- 
ments are not the same for all Ti/0 atoms in case of 
[110] strain. Two inequivalent Ti (Ti(l), Ti(2)) and O 
(0(1), 0(2)) sub-lattices exist under such strain, as the 
equatorial Ti(2)-0(2) bonds are not modified while the 
imposed strain enlarges the apical Ti(2)-0 bonds (the 
opposite holds for the Ti(l)-0 distances), see Fig. [l](b). 
As a result, the short-range repulsion and thus the ampli- 
tudes of the atomic displacements differ for both atomic 
subclasses. Specifically, the displacements along [001] are 
larger for Ti(l)/0(1) than for Ti(2)/0(2) (and the other 
way round for the [110] displacements), cf. Refs. ITTland 



As already mentioned previously, the energy difference 
between the paraelectric and the ferroelectric phase is 
often discussed in terms of a delicate balance between 
short-range repulsive forces and the long-range dipolar 
interaction which favor the ferroelectric state.^^ It has 
been further shown for perovskites like BaTiOa that 
the dipolar interaction can counterbalance the short- 
range repulsion if an anomalous hybridization between 
the displaced ions lowers the energy during a ferroelectric 
transition.'^ The local chemical environment in Ti02 is 
quite similar to the perovskite structure, with Ti^+ ions 
that are octahedrally coordinated with 0^~ ions, and 
mostly ionic bonds with some covalent character. It is 
thus very likely that similar dynamic hybridization effects 
as in perovskites can also stabilize the ferroelectric state 
in strained rutile. Indeed, a Mulliken population analy- 
sis showed an increased Ti^^-Op hybridization within the 
ferroelectric phase.^^ 

To obtain a detailed understanding of this dynamic 
charge transfer, we make use of the dynamic Born effec- 
tive charge. 



dPi 
\e\ dr^. 



(3) 



which describes the change in polarization P in direction 
i for a shift of atom hi in direction j. Here, Q is the 
unit cell volume and |e| the electronic charge. Due to 
the charge neutrality condition, 2 • Z^- = —A-Zq, we 
can restrict the following discussion to Z^-. The effec- 
tive charge tensor of Ti in rutile is diagonal in a [110], 
[001], [110] reference system, see Fig.[l](b), with the prm- 
cipal values ^*oqi], and Z^. Here, Z|* is the dynamic 
charge along the apical bond, which is pointing along 
[110] respectively [110] for every second Ti-atom, and 
is the charge along the third perpendicular direction in 
each case. 

In Tab. |l] our values calculated with VASP and the 
values obtained form the shifts of the center of gravities 
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TTfictl nfiCT TB*2i TSjS 

(a) (b) ^^^^^ ^^^it 

FIG. 3. (Color online) Maximally localized Wannier functions corresponding to (a) sp^g (b) splq^^p^ (c) spap and (d) 
p^-orbitals. 



of the maximally localized Wannier functions (MLWFs) 
based on PWscf calculations are compared to literature. 
Although the VASP results are slightly larger than all 
other values, due to the use of soft potentials, the quali- 
tative trends are reproduced within both approaches. 

It can be seen that in particular the [001] -component 
of the Born tensor, ^[qqi]' exceptionally large and 
of the same size as for ferroelectric materials (BaTiOs: 
Z^Q-i^l=7.3).'^ Furthermore, the Born charge along the 
apical bond, is also quite large and may stabilize a 
ferroelectric transition in [110] /[1 10] direction, whereas 
is only moderately enhanced compared to the nomi- 
nal ionic charge of +4. 

As a first step towards the understanding of dipolar 
interactions at the rutile surface it is essential to under- 
stand in more detail the underlying mechanism for the 
exceptionally large Born effective charge of Ti'^^ in bulk 
Ti02, and how it is affected by different lattice modifi- 
cations. These aspects will therefore be discussed in the 
remainder of this section. 

Based on the constructed MLWFs we decompose Z^- 
according to (see, e.g., Ref . [42]) : 

z* = z:^,,+z:^+z:^i , (4) 

where Z^^^^ = +12 is the charge (in units of |e|) of the 
nucleus screened by the inner electrons, Z*^ the contri- 
bution of the Ti 35 and 3p semi-core states, and Z*^^ is 
the contribution of the valence electrons. We find that 
the semi-core contribution to Z*^ [ooi]) i^sc ||) ^i^^^s only 
by —0.23 (—0.12) from the nominal value of —8. In con- 
trast, large anomalous values appear for Z*^^, which can 
be further decomposed into contributions of individual 
MLWFs, see Tab. [V[ Thereby, the 16 valence bands of 
Ti02 with dominant 02s and 02p atomic character are 
represented by 16 MLWFs per unit cell, 4 orbitals per 
O ion. Starting from projections onto three atomic sp'^ 
hybrid orbitals, centered at each of the four different oxy- 
gen atoms and oriented towards the three surrounding Ti 
atoms, and one p orbital oriented perpendicular to the 
corresponding O-Tia triangles (see Fig. [TJa)), we obtain 
three a and one 7r-type orbital per oxygen as depicted 
in Fig. [3j We note that depending on the initial projec- 



tion used for the spread minimization, different localized 
Wannier orbitals can be obtained, corresponding to vari- 
ous local minima of the total quadratic spread functional. 

In agreement with the formal ionic Ti^+-0^~ configu- 
ration, the centers of gravity of each MLWF are close to 
the central 0-atoms. Nevertheless, the mixed covalent- 
ionic character of the Ti-0 bonds can clearly be seen 
from the d-character "tails" of the MLWFs at the Ti po- 
sitions (see Fig.|3|. This is consistent with the integrated 
projected density of states shown in Fig.js) which, for ex- 
ample, results in an occupation of about 0.6 electrons for 
the Ti dxy atomic orbital.*^ Furthermore, the center of 
gravity of the a-type MLWF along the shorter equatorial 
bond divides the bond with a ratio of 22:78, while the 
longer apical bond is divided with a ratio 20:80. 

If the central Ti-atom is displaced along [001] ([110]), 
the Ti-Oeg (Ti-Oap) distances are modified, and thus the 
covalent character of the corresponding bonds change. 
This leads to a shift of the center of gravity of the ML- 
WFs, which causes large and positive anomalous contri- 
butions to the Born effective charges, see Tab. |V| (we 
note that in the fully ionic case the valence contribu- 
tion to the Ti Born effective charge would be identical 
to zero). Thus, charge is dynamically transferred in and 
out of the Ti-d orbitals. Similar to the case of the per- 
ovskite ferroelectric BaTiOs, the 7r-type orbitals {p±) 
show the largest anomalous contributions to Z*, under- 
lining the large similarity between Ti02 and the related 
ferroelectric perovskites. 

In addition to these positive contributions, the ML- 
WFs located at Oap (Ogg), exhibit small negative con- 
tributions to Zj^QQ-^i (^1*)* "^^^ corresponding Ti-0 bond 
lengths do not change to first order in the displacements 
and thus no strong change in hybridization occurs. In- 
stead, the charge contained in the "tails" of the corre- 
sponding MLWFs move with the Ti atom, leading to 
small negative contributions to Z*. This effect is espe- 
cially pronounced for the equatorial 7r-type orbital for a 
displacement of the Ti atom along [110], thereby slightly 
reducing Z|*. 

We note that a similar Wannier decomposition for the 
Born effective charge of the oxygen anion in rutile Ti02 
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4. (Color online) (a)-(b) Change in the Ti Born charge 
(a) and Z*| (b) due to a modification of the equatorial 
Ti-0 bond length by uniform lattice expansion (red, squares), 
modification of u (blue, diamonds), [001]-strain (black, trian- 
gles), and [110] strain. For the case of [110] strain the given 
Ti-Oeq distance corresponds to Ti(2) and the Born charges 
are depicted for both Ti(l) (cyan, stars) and Ti(2) (green, 
pluses) . 



has been presented in Ref. [43| In agreement with our 
work, Cangiani et al found the largest anomalous con- 
tributions to Zq for the p^-type orbitals. However, even 
though the total charges of Ti and O are related by the 
charge neutrality condition, the relation between indi- 
vidual orbital contributions is less straightforward, thus 
preventing a more detailed comparison. 

One important factor for the modification of the dipo- 
lar interaction at surfaces are atomic relaxation and 
structural distortions at the surface. In the following we 
therefore consider different lattice modifications within 
the bulk system and investigate the resulting effect on 
the Born effective charges. Thereby we are considering 
similar lattice modifications as in the previous section, 
i.e., a uniform lattice expansion, modification of the ra- 
tio of the Ti-0 distances by and imposed strains in 
the [001] and [110] directions, in each case without fur- 
ther relaxation of the other parameters. The results are 






+ 


dyz=dxz 




dz2-r2 








Px=Py 




Pz 





FIG. 5. (Color online) Orbitally-resolved changes in the in- 
tegrated valence density of Ti states due to a modification of 
the Ti-0 bond lengths. Solid lines: Modification of the inter- 
nal parameter u. Dashed lines: Uniform lattice expansion, x: 
[001], y: [010], z: [001]. 



depicted in Fig. [4j 

A uniform lattice expansion leads to a small increase 
of both 'Zj^qqi^ and while at the same time the static 
hybridization of the atoms is reduced, see Fig. [5] For 
all other distortions the ratio of the Ti-0 bond lengths 
is modified in addition to eventual changes in the vol- 
ume. A variation of the internal parameter u in turn con- 
serves the volume, while increasing the equatorial bond 
distance and decreasing the apical bond distance. As 
seen in Fig. |5] this also changes the static hybridization. 
On one hand, the occupation of the Ti-d and p states 
within the x-y plane increases with the largest gain for 
the dxy orbital, which is oriented along the apical bond. 
On the other hand, the occupation of orbitals with lobes 
in [001] direction is reduced. Thus, the Ti-Ogg [Ti-O^p] 
hybridization decreases [increases], as do the weights of 
the corresponding tails in the MLWFs. 

This change in the static hybridization seems to reduce 
the dynamic change in hybridization between the Ti and 
Oeq orbitals especially for the p± orbital. As a conse- 
quence, both the positive and negative contributions of 
this orbital to 'Zj^qqi^ and Z*|, respectively, are reduced. 
Similarly, the reduced coupling of the sp^^p orbital with 
the equatorial Ti neighbor may explain the slight reduc- 
tion of its contribution to 'Zj^qqi^ In contrast, the hy- 
bridization between the Ti atom and its apical neighbor 
increases for increased as do most of the negative and 
positive contributions of these orbitals to ZJqq-^| and Z*|. 



In summary, a systematic reduction [increase] of the ef- 
fective charge along [001] ([110]) occurs, see Fig. |4] For 
example, Z*^^ ^qq-^^ i^vai ||] decreases [increases] by 8 % 
[24 %] within the PWscf calculation for a u modification 
of +3 %. Thereby, the main reduction in Z*^^ ^qq-^^ is due 
to a reduced anomalous contribution of the p± orbital at 
Oeq, whereas the increase of Z*^^ n is due to slightly en- 
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larged contributions of pj_ and sp^^ at Oap and a strong 
reduction of the negative contribution of p± at Oeg. 

One may wonder, whether the absolute Ti-Oap dis- 
tance instead of the Ti-0 bond ratio plays the major role 
for this modification of the static and dynamic hybridiza- 
tion. However, if the Ti-Ogg distances are modified for a 
fixed Ti-Oap distance by a combination of isotropic ex- 
pansion and a variation of the modification of ^[qq^] 
with the ratio of the Ti-0 distances is equal to a variation 
of u. 

As the change of the Ti-0 bond length ratio is smaller 
in the case of [001] strain, the corresponding reduction 
[increase] of ^[qqi] l^vai ||] ^^^^ pronounced. Similar to 
the modification of Z*|, the third principle value of the 
Born tensor, slightly increases with a reduced ratio 
of the Ti-0 distances. However, for small variations of 
the Ti-0 distances, as they may appear at the surface, 
the absolute value of is not significantly enhanced 
compared to the formal ionic charge, and thus no con- 
siderable infiuence on the ferroelectric trends has to be 
expected from this modification. 

For [110] strain, the ratio of Ti-0 bonds correspond- 
ing to Ti(2) is modified analogously to the case of [001] 
strain, and indeed the corresponding decrease of ^[qoi] 
is identical to [001] strain for small strain values, see 
Fig. Qa). In contrast, the apical bond length increases 
for Ti(l) while the equatorial distance stays constant, 
and thus Z^^^-^^ increases for this Ti atom. Furthermore, 
while Z|* is nearly independent of the equatorial bond 
distance for Ti(2) under [1 10] strain, Z|* decreases dras- 
tically for Ti(l) with the corresponding expansion of the 
apical bond. 

In summary, the largest individual contributions to the 
anomalous Born effective charges in rutile stem from the 
large dynamic charge transfer along Ti-0 tt bonds. Ad- 
ditionally, the static Ti-0 hybridization, which depends 
strongly on the ratio of both Ti-0 distances, leads to 
considerable modifications of the effective charges. The 
resulting modification of the dipolar interaction has a 
considerable influence on the ferroelectric properties of 



the rutile surface, as we will see in Sec. IV. 3 



Furthermore, we find large effective charges which are 
stable against small lattice distortions, as they may ap- 
pear at the rutile surface. 



IV. (IIO)-SURFACE 

For the investigation of the ferroelectric trends of the 
rutile surface, we use two different approaches. After 
summarizing the structural relaxations at the surface we 
first briefiy discuss the ferroelectric trends for a free film 



of 7 monolayers (ML), whereas in Sec. IV. 3 and IV. 4 



model the clamping of the film due to an idealized sub- 
strate by fixing two bottom layers to paraelectric, respec- 
tively ferroelectric positions of the bulk material under 
different strain conditions (see Tab. IV). 



(110) 




(110) 



Layer 3 



FIG. 6. (Color online) Atomic structure of the topmost layers 
of the rutile (110) surface. Light cyan: Ti, red: O. Black 
arrows: Inter layer distance oscillation. 



As can be seen in Fig.[6j the (110) surface consists of 
Ti-Op layers with bridging oxygen (O5) atoms between 
these layers. The surface Ti-atoms are alternately five- 
fold coordinated (Tisc) and six-fold coordinated (Tiec). 
We note that we use the nomenclature Ti^c (Tisc) not 
only for the atoms in the first surface layer, but also for 
the corresponding Ti atoms below. 



IV. 1. Surface relaxation 

The calculated atomic structure of the paraelectric 
films is very similar for free and clamped films and thus 
we restrict the discussion to the structural relaxation of 
the clamped film. The corresponding results are sum- 
marized in Table |VIl Even though the rutile (110) sur- 



face is non-polar and stable without reconstruction, the 
atomic structure is significantly modified at the surface. 
It has been discussed in literature that the atomic struc- 
ture converges only slowly with film thickness.'^'^'^ How- 
ever, we yield modifications of the obtained interatomic 
distances smaller than 0.01 A if the increase the film 
thickness to 9, only. Additionally, our obtained struc- 
tural relaxation at the su rface agrees well with previous 
results in literature) ^ * ^ * -'^^ * 

Most notably, the outermost Ti^c and O5 atoms in the 
first surface layer relax inwards and reduce their distance 
towards the underlying apical and equatorial neighbors 
by 9 %, respectively 6.1 % (6.6 %), in our calculations (in 
Ref. ^). In contrast, the next deeper Ti-0 bonds along 
the surface normal are enlarged by 3.6 %, respectively 
7.1 %, for the equatorial Tisc-O^ and Tiec-O^ pairs. This 
alternating modification of the Ti-0 bonds reaches sev- 
eral A deep into the surface, cf. Refs. 7 and iTOl 

Apart from this oscillation of the Ti-0 distances, a 
small surface buckling appears, as Tisc-atoms and bridg- 
ing 0-atoms (O5) relax inwards whereas Tiec and O 
atoms within Ti-0 planes (Op) relax outwards, see Fig. [6] 
It has been debated in literatur^^ whether an alike polar 
surface relaxation of SrTiOs can be interpreted as ferro- 
electric relaxation along the surface normal. However, 
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TABLE VL Ti-0 bond lengths (A) within the three outermost layers of the rutile (110) surface within the paraelectric state 
(see Fig.[6| with and without strain. The bulk distances refer to the optimized paraelectric structure with corresponding strain. 
Left columns: Tiec: sixfold coordinated Ti-atoms in the surface and Ti-rows below. Right columns: Tisc: fivefold coordinated 
Ti-atoms in the surface and Ti-rows below. A single entry corresponds to Ti-0 neighbors in plane while each Ti atom possesses 
different distances towards its inner/outer Ob neighbors. 
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this relaxation does not break the symmetry of the ideal 
surface and thus cannot be interpreted as a ferroelectric 
transition, nor exist two switchable states.*^ 

In summary, the Ti-0 bonds along the surface normal 
are considerably modified, and show a large anisotropy 
for each Ti-atom, whereas the bond lengths perpendic- 
ular to the surface normal are very close to their bulk 



values, see Tab. VI 



IV. 2. Free surface 

Due to the strong symmetry lowering at the surface, 
the ferroelectric phonon modes of the bulk are no longer 

dis- 
This 



eigenmodes of the system. Most notably, the Ti -Oe q 
tances vary between 1.83 A and 2.11 A, see Tab. VI 



means that the short-range repulsion during a ferroelec- 
tric shift is locally very different, and one may expect 
different amplitudes of polar displacements in each layer. 
As a result of this, a polar displacement of the whole slab 
along the bulk A2n mode (i.e., with the same relative dis- 
placement between Ti and O in all layers) is energetically 
very unfavorable. The stiffness of such a displacement in- 
creases by a factor of ten in comparison to bulk. 

However, if all atomic positions are optimized after 
such a static displacement along the bulk A2n mode, 
small polar displacements persist within the film. For 
instance, the total energy of the paraelectric system is 
lowered by 3 meV/cell if mainly the Tigc atom in layer 3 
is displaced by 0.07 A along [001] against its surrounding 
O octahedra, resulting in a dipole moment of 0.8 |e|-A 
(determined with Eq. [s]). For this Ti atom, the enlarged 
short-range repulsion due to the shrink of two equatorial 
bonds of 2 % can be compensated, as the apical bond 
is nearly conserved and the two other equatorial bonds 
increase by about 3 %. A ferroelectric configuration be- 
comes even more favorable if small polar displacements 
in [110] direction are superimposed to the displacement 



along [001], and we obtain an energy gain of 5 meV/atom 
relative to the paraelectric surface. 

Although these small energy well depths will most 
likely be overcome by thermal or even quantum fluctua- 
tions, the presence of such polar surface states confirms 
the strong ferroelectric trend at the rutile surface, which 
is comparable to the SrTiOs surface.'^ In addition, these 
polar displacement patterns demonstrate the large modi- 
fication of the polar phonon mode at the surface. Thus, in 
order to investigate the ferroelectric trends of the free sur- 
face using the frozen phonon approach (similar to what 
has been done for the bulk material, see Sec. |IIL2 ), it 
would be necessary to calculate mode eigenvectors for 
the strained films. However, since the determination of 
surface phonon modes is computational very demanding, 
we will follow a slightly different approach in order to 
investigate the ferroelectric trends of the surface in the 
remainder of this paper. 



IV. 3. Clamped surfaces 

We use the atomic structure of the bulk material within 
the paraelectric or ferroelectric phases, see Tab. [Hj as 
starting point and fix two bottom layers to these bulk 
positions. This setup corresponds to a clamping of the 
bottom layers to an underlying substrate. By using this 
procedure, the amplitude of the displacements within the 
upper layers can account for the local atomic arrange- 
ment of the surface, and thus will approach the real eigen- 
state of the system. 

To quantify the local polar displacements, we use the 
total dipole moment pi per Ti204 layer i: 



(5) 



where Z*-^ is the Born charge of atom j along r, and 
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FIG. 7. (Color online) Layer resolved Ti density of states: 1st layer (yellow, light, dotted); 2nd layer (red, darker, dotted); 3rd 
layer (blue, dark, solid line); 4tli layer (cyan, light, solid line) and bulk value as thin black lines, (a) The numbering starts at 
the fixed bottom of the film (b) The numbering of layers starts at the relaxed surface. 



Afj is the shift with respect to the paraelectric reference 
configuration of each atom j within layer i. 

As mentioned previously, the surface properties of 
Ti02 converge only slowly with film thickness. This is 
especially important for electronic properties such as the 
band gap.^^^^ Furthermore, the clamping of one surface 
of the slab could induce artificial electronic surface states 
which could then modify the ferroelectric trends. To 
address this issue, Fig. [7| shows the layer-resolved total 
density of states (DOS) corresponding to the Ti atoms. 
Indeed, surface states appear within the band gap at 
the fixed bottom surface, see Fig. [7|^a). However, these 
states decay rapidly with increasing distance from the 
bottom layer and do not significantly influence the elec- 
tronic structure of the free top surface. In addition, the 
DOS obtained here is in qualitative agreement with the 
electronic structure obtained in previous investigations.^ 
The largest modification of the electronic structure in 
comparison to the bulk material is a slight narrowing of 
the band gap, due to a small shift of the unoccupied d- 
states of the Ti^c atoms. 

As discussed by Bredow et al.^ the hybridization be- 
tween Ti and O states which are oriented along the sur- 
face normal, oscillates with the surface distance. These 
surface-induced modifications of the hybridization are 
correlated with the reduced (enlarged) Ti-0 distances at 
the surface, see Sec. |IV.1| As we showed in Sec. |III.3[ 
small changes in the Ti-0 arrangement and hybridiza- 
tion can have considerable influence on the ferroelectric 
characteristics. Because of this, it is essential to have 
a detailed understanding of the character of the Ti-0 
bonding in the vicinity of the surface in order to make 
predictions on the ferroelectric trends. For this purpose, 
we decompose the valence states of the surface slab into 
MLWFs, see Fig.jS] We note that in order to obtain qual- 
itative trends at a reasonable computational effort, this 
decomposition has been performed for a free film with a 
thickness of only three monolayers. 

It can be seen from Fig.[8];a)-(b) that the MLWFs with 
sp'l and p± character which belong to the topmost O5 



atom are considerably modified in comparison to bulk. 
For both orbitals, the Ti-atom in the apical position is 
missing at the surface, and the centers of gravity shift 
towards the Ti^c atom below, see subfigure (c). Thus, 
the covalent character of this bond increases. Similarly, 
the hybridization of the undercoordinated Ti^c surface 
atom and the O5 atom below is enlarged. Overall, an 
alternating increase and decrease of the Ti-0 hybridiza- 
tion appears along the surface normal, and the amplitude 
of this oscillation decreases with the surface distance, 
fully consistent with the changes in the corresponding 
bond lengths. Thereby, the strongest shifts correspond 
to the centers of gravity of the p^-type MLWFs, which, 
as shown in Sec. III.3[ also exhibit the largest individual 
contributions to the bulk Born effective charges. Further- 



more, we also found in Sec. III. 3 that a modification of 



this p^ orbital also leads to a considerable modification 
of the dynamic charges, and thus an oscillating modifica- 
tion of the dynamic charges can be expected at the rutile 
surface. 

In order to verify this conclusion, we now discuss the 
modifications of the Born effective charges at the (110) 
surface in more detail. Fig. |9] shows the Born charges 
along [001] and [1 10] for each class of atoms. ^ We only 
consider cases with polarization within the surface plane, 
in order to avoid complications due to the strong depo- 
larizing fields perpendicular to the surface. The Born 
charges are calculated both for the structurally relaxed 
as well as for the unrelaxed surface configuration. 

For a detailed discussion of dipolar interactions along 
[110], one has to keep in mind that Z*| and are al- 
ternating in this direction. While all apical bonds are 
aligned along [110] for O^ (Z*j-^q|=Z^), they are aligned 
within the surface for the Op atoms (Z^j-^q|=Z*|). For the 
Ti- atoms, Z*| and Z^ are alternating in each layer, see 
Fig.pd-f). 

It is apparent from Fig. [9] that even if atomic relax- 
ations are neglected, the dynamic charges are strongly 
modified within the topmost layers. In particular, ZJqq-^| 
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(c) 



FIG. 8. (Color online) (a)-(b) MLWFs of the topmost Ob atom, (a) sp^p orbital (b) pJ_ orbital, (c) Centers of gravity of 
the MLWFs close to the surface. Light, yellow spheres: spe^; dark, blue spheres: sp^^; black spheres: p^-orbital. The bulk 
distribution of the MLWF centers are given as reference within the small boxes. Arrows mark the overall shift of charge at the 
surface relative to the bulk structure. 
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and Z^j-^Q| of Op [(e)] atoms are slightly 
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Ti-atoms [(f)] 
reduced. 

If we account for atomic relaxations, the alternating 
strong and weak Ti-0 bonds appear along the surface 
normal, as discussed above. Especially, the modification 
of the Ti-Oeg bonds has a considerable effect on the Born 
charges along [001]. Simultaneous with the modification 
of the equatorial Ti-O^ bonds, Z^qq-^^ of these O5 atoms 
oscillates, i.e. the dynamic charge transfer increases (de- 
creases), if the corresponding Ti-Oeg bonds are strength- 
ened (weakened), while the amplitude of this oscillation 
decreases with increasing distance from the surface (as 
does the oscillation of the bond strength). This is fully 
consistent with the reduction (increase) of Z^qq-^^ for a re- 
duced (enlarged) Ti-Ogg static hybridization of the bulk 
material, see Sec. |III.3 



FIG. 9. (Color online) Layer-resolved Born charge at the ru- 
tile (110) surface, (a)-(c) ZJqq^^] (d)-(e) Ta^iiqy Bulk values of 
Zfooi]5 Z*| and Z^, (the mean value of ZJ^j^q^) are marked by 
solid (dotted) lines, (a) /(d): Ob (b)/(e): Op, (c)/(f): Ti; 
Tiec atoms and the corresponding Oeq atoms (red, diamonds); 
Tisc and the corresponding Oeg atoms (blue, cross); mean 
value (cyan, triangles); mean value of the unrelaxed cell 
(black, squares). Lines are a guide to the eye only. 



of the outermost Tiec [(a)] and O5 atoms [(c)] is signif- 
icantly enlarged in comparison to bulk, while the cor- 
responding charges within the topmost Tisc [(a)] atom 
and its equatorial Op [(b)] neighbors is slightly reduced. 
Also, the dynamic charges along [110] are reduced within 
the topmost surface layer, as the mean value Z^^-^^^ of the 



We note that the charge neutrality condition does not 
hold for each layer under the reduced symmetry at the 
surface, and as a result about 0.2 electrons are trans- 
ferred from the dynamic charges of layer 1 to the ones of 
layer 2. In particular, the Born charges of the uppermost 
Tisc atom and the Tigc atom in layer 2 are not balanced 
by their equatorial neighbors. For these Ti atoms, the 
Ti-Oap bonds along the surface normal are reduced to 
1.83 A, respectively 1.87 A and are thus shorter than the 
bulk Ti-Ogg bonds. Because of this, the Born charge in 
layer 1 is slightly reduced. Besides the direct surface, 
the mean Born charges per layer of Ti, 0^ and O5 atoms 
oscillate with a slight increase in odd layers. 



The mean values of Z; 



110] 



also oscillate with the sur- 



face distance for Ti and Op atoms, see Fig. [9]^e)-(f), but, 
the Born charge is enlarged within even layers opposite 
This behavior is consistent with the opposite 
and Z*| for different modification of the 



^fooi]- 
trends of Z^qq-^^ 
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[110]-strain 
(b) 



4 6 
Layer 



FIG. 10. (Color online) Layer-resolved dipole moments under 
tensile uniaxial strain, (a) |p[ooi]| (b) |P[iio]|- Solid (dotted) 
lines correspond to strain parallel (perpendicular) to the po- 
larization. 2 % compressive strain (yellow, stars); no strain 
(black, crosses); 2 % tensile strain (red, squares); 5 % tensile 
strain (blue, diamonds). In all cases polar distortions per- 
pendicular to the polarization direction of the substrate are 
prevented by the imposed symmetry. 



Ti-0 bond ratio and hybridization in the bulk. of the 
O5 atoms is nearly constant while of Tisc C^hc) in- 
crease (decrease) marginally. As the Ti-0 bonds within 
the surface planes are modified to a smaller extent, the 
oscillations of Z*| of 0^ atoms are only very minor. Be- 
sides this, Z*| of Tiec C^hc) is reduced (enlarged) in all 
layers and a dynamic charge transfer of about 0.2 |e| ap- 
pears from layers with odd numbers to layers with even 
number. 

Overall, besides the discussed modifications of the Born 
effective charges, the large anomalous dynamic charges 
along [001] and [110] found in the bulk system are re- 
tained at the (110) surface of rutile. Therefore, a similar 
tendency towards a ferroelectric distortions as in the bulk 
material can also be expected at the surface. 



IV. 4. Strained surface 

In the following, we present our results on polar atomic 
distortions of the clamped surface. For the unstrained 
film we tested the stability of polar distortions along 
[001] by a displacement of 0.1 A along the polar 



bulk mode of the 2 clamped bottom layers, cf. Sec. III. 2 



As the upper layers are fully relaxed, this setup corre- 
sponds to the clamping of the Ti02 film to a ferroelec- 
tric substrate. Although, small polar displacements are 
stable with respect to the paraelectric configuration in 
layer 3-5, the amplitude of these polar displacements de- 
cays with the distance towards the polar layers and thus 
no stable polarization is obtained directly at the surface, 
see Fig. IlOfa). This decline of the polarization shows 




FIG. 11. (Color online) Layer-resolved dipole moments under 
tensile uniaxial strain for an optimization of the atomic struc- 
ture without imposed symmetry. The substrate layers have 
been fixed to the corresponding polar states along (001). Solid 
lines: |p[ooi]|5 dashed lines: |p[iio]|- 2 % strain (red, squares); 
5 % strain (blue, diamonds); (a) [001] strain; (b) [110] strain. 



that the ferroelectric state is not stable within the un- 
strained film although a large polarizability is present. 
Furthermore, the decay of the polar distortions is even 
more pronounced if compressive [001] strain is imposed. 
This indicates that a ferroelectric distortion might be sta- 
bilized by a systematic expansion of the Ti-0 distances, 
cf. discussion of the bulk case in Sec |III.2| 

Indeed, polar distortions along [001] or [110] are stable 
within the Ti02 surface for 2 % uniaxial tensile strain, 
see Fig. [lO] Additionally, the short-range repulsion is 
further reduced for an increasing tensile strain and thus 
the local dipole moments increase systematically if we 
impose an uniaxial expansion of 5 %. 

Table |VII| lists the energy gain for the different polar 
displacements in comparison to the strained paraelectric 
films. In agreement with the bulk system, both polar 
states along [001] and [110] are energetically more fa- 
vorable than the paraelectric state and thus two stable 
ferroelectric states exist. Obviously, the polar distortions 
parallel to the imposed strain are energetically most sta- 
ble. But, in agreement to the strained bulk, the para- 
electric state is also unstable against a polar distortion 
in the direction perpendicular to the strain, see Tab. |IV| 
Only for 2 % [001] strain the paraelectric state is more fa- 
vorable than the ferroelectric state perpendicular to the 
strain direction. 

For the obtained dipole patterns in Fig. [Toj we con- 
strained the atomic optimization and enforced the local 
dipoles to point parallel to the fixed bottom layer. In 
contrast, if all atomic degrees of freedom are optimized 
without such constraint, local dipoles both along and per- 
pendicular to the polarization of the fixed bottom layer 
For the polar states along [001], the 



form, see Fig. 11 



system obviously gains energy, if superimposed polar dis- 
tortions along [110] are formed. Tab VII The energies for 
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TABLE VII. Energy gain relative to the corresponding para- 
electric case (meV/2 f.u. surface cell) of the different polar 
surface configurations, see text. The bottom layers are fixed 
to the polar state along [001] ([110]) in the first (second) row. 
Only in-plane displacements parallel to the bottom polariza- 
tion are allowed. (The upper layers are fully relaxed resulting 
in local dipole moments along [001] and [110]). 

[001]-Strain [nO]-Strain 
Bottom layers 2 % 5 % 2 % 5 % 

[001] polarized 18 (25) 262 (275) 3 (48) 111 (518) 

[110] polarized -1 48 (128) 28 525 



polar distortions in both directions are given in brack- 
ets. The modification of the local dipoles is even more 
pronounced under tensile strain along [110]. Polar distor- 
tions along [001] can be stabilized under such strain if the 
symmetry of the system is fixed, which could possibly be 
achieved by an external electric field. Despite this, polar 
distortions parallel to the strain direction are much more 
favorable and are even stable on top of a ferroelectric 
substrate polarized along [001], see Fig. 11 'b). 

We now discuss the observed layer-dependent varia- 
tions of the local dipoles. In all cases, the local dipole 
moments in [001] ([110]) direction are enlarged in lay- 
ers with an odd (even) number. Besides this, the local 
polarization in [001] direction slightly decreases towards 
the surface and is quenched within the first first layer for 
[110] strain or 2 % expansive [001] strain. For the local 
polarization along [110] no systematic decay of the po- 
larization towards the surface is obtained, although, the 
polarization in the first layer is considerably reduced for 
all kinds of strain. 

The different Ti-0 distances and thus the local vari- 
ations of the short-range repulsion have a large influ- 
ence on the obtained dipole patterns, see discussions in 
Sec. [no and Sec. [TO The Ti-0 distances at the sur- 
face and thus the short-range repulsion oscillates with 
the distance form the surface. Additionally, the tensile 
strain further modifies the spatial variation of the Ti-0 
distances. First of all, the Ti-Ogg distances within the 



Ti-Op planes increase systematically with both kinds of 



tensile strain, see Tab. VI In addition, the apical dis- 
tances in-plane increase with [110] strain whilst they are 
not modified by strain along [001]. Because of this, the 
short-range repulsion which builds up during a polar shift 
between Ti atoms and their Op neighbors in [001] and 
[110] direction is reduced considerably and thus polar 
distortions are stabilized by tensile strain. 

As the tensile in-plane strain enlarges the atomic vol- 
ume, a shrink of the atomic distances along the surface 
normal is very likely. Indeed, the Ti-Ti interlayer dis- 
tances decrease for tensile [001] strain. In contrast to 
this, no unique reduction of the interlayer distances can 
be found for tensile [110] strain. 

Additionally, for the Ti-0 bonds oriented perpendic- 
ular to the surface, the discussed strengthening (weak- 



ening) of the bond has a significant impact on the 
atomic relaxation under imposed strain. On one hand, 
the strengthened bonds are rather insensitive to tensile 
strain. In particular, the distance between the Tisc atom 
in layer 1 and the outer O5 atom is not modified by 2 % 
[001] strain or [110] strain and shows an expansion of only 
1 % under 5 % [001] strain. As a consequence, a polar 
displacement of this Tisc atom relative to the surround- 
ing atoms along [001] is not stable, even at 5 % expansive 
strain in either direction. On the other hand, the length 
of the less covalent bonds increases (decreases) system- 
atically under expansive [001] ([110]) strain. Because of 
this, the asymmetry of the Ti-0 bond lengths along the 
surface normal increases (decreases) slightly with the im- 
posed [001] ([110]) strain for each Ti atom. 



Besides the modified short-range repulsion, the dy- 
namic charges at the rutile surface oscillate with the dis- 
tance from the surface, see Fig. |9] Consequently, the 
local dipole moments show the same trends with maxi- 
mal (minimal) local dipole moments along [001] in odd 
(even) layers and the other way around for polarization 
along [110]. Here, the magnitude of Z* has a direct in- 
fluence on the local dipole moments, see Eq. [5] Also, the 
local variation of the dipolar interaction leads to a mod- 
ulation of the atomic displacements. 
Thus, the considerable difference between Z*| and 
leads to different amplitudes of the [1 10] displacements 
of Tisc and Ti^c atoms in each layer. For example, 
the apical bond of the Tiec atoms in the first layer is 
aligned along [110] and thus the displacement of this 
atom against the surrounding oxygen octahedra is a fac- 
tor of two larger than those of the Ti^c atom within the 
same layer. Analogously, the displacement of the Tisc 
atom in each even layer is a factor of two larger than the 
displacements of the corresponding Tiec atom. Addition- 
ally, Z|| is considerably larger within the even layers, as 
are the amplitudes of the polar distortions. Similarly, the 
reduction of P[ooi] the layers with an even number is 
due to a reduced displacement of the Tiqc atoms within 
these layers, which is correlated with the reduced dipolar 
interaction along [001] for these atoms. 



Both Zfoo,] 



sile strain, see Tab. 



and Z^^Q^are slightly modified under ten- 



VIII 



For example, Z|qq-,^| of the O, 



atoms increases (decreases) for the strengthened (weak- 
ened) bonds, while the Born charge of the Op atoms is 
insensitive towards such strain. As a consequence, the os- 
cillation of the dipolar interaction is slightly enlarged in 
comparison to the equilibrium lattice constant. For [110] 
strain 'Zj^qqi-^ of all Op and Ti atoms is enlarged and thus 
the magnitude of the Z|qq-^| oscillation is reduced. Simul- 
taneously, the layer- wise oscillation of P[ooi] is smaller for 
[110] strain. 
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TABLE VIII. Born charges (|e|) within the paraelectric phase of the rutile (110) surface calculated with VASP for different 
values of uniaxial [001] and [110] strain. For the Tiec/Tisc atoms and the Ti atoms below, respectively the Op and Ob atoms, 
the Born charge for the first free surface layers are given, see Fig.pj Upper part ZJqq^^]; Lower part: ZJ'j^^qj. 







Tiec 






Tisc 






Op 






0, 




Strain 


1 


2 


3 


1 


2 


3 


1 


2 


3 


1 


2 


3 


No strain 


8.4 


7.3 


8.3 


7.3 


7.6 


8.2 


-3.8 


-3.7 


-4.1 


-5.6/-2.8 


-2.8/-4.5 


-4.6/-3.7 


[001] +2% 


8.6 


7.2 


8.3 


7.3 


7.6 


8.2 


-3.8 


-3.7 


-4.1 


-5.8/-2.7 


-2.7/-4.6 


-4.7/-3.6 


[001]+5% 


8.7 


7.1 


8.4 


7.2 


7.6 


8.1 


-3.8 


-3.6 


-4.1 


-6.1/-2.6 


-2.5/-4.7 


-4.8/-3.5 


[iio]+2% 


8.6 


7.3 


8.4 


7.4 


7.7 


8.3 


-3.8 


-3.7 


-4.2 


-5.7/-2.9 


-2.9/-4.5 


-4.7/-3.8 


[iio]+5% 


8.9 


7.5 


8.7 


7.6 


8.0 


8.4 


-3.9 


-3.8 


-4.2 


-5.7/-3.1 


-3.1/-4.6 


-4.7/-4.0 


No strain 


6.3 


5.7 


7.2 


4.5 


7.8 


5.1 


-4.3 


-5.5 


-4.9 


-1.2/-1.3 


-1.1/-1.2 


-1.4/-1.4 


[001]+2% 


6.4 


5.7 


7.4 


4.6 


7.9 


5.1 


-4.4 


-5.6 


-5.0 


-1.1/-1.3 


-1.1/-1.2 


-1.4/-1.4 


[001]+5% 


6.7 


5.9 


7.6 


4.7 


8.0 


5.3 


-4.6 


-5.6 


-5.2 


-1.2/-1.3 


-1.1/-1.2 


-1.4/-1.4 


[iio]+2% 


6.1 


5.7 


7.0 


4.4 


7.8 


4.9 


-4.1 


-5.5 


-4.8 


-1.1/-1.3 


-1.1/-1.2 


-1.4/-1.4 


[iio]+5% 


5.9 


5.7 


6.7 


4.2 


7.6 


4.8 


-4.0 


-5.3 


-4.6 


-1.2/-1.3 


-1.1/-1.2 


-1.4/-1.3 



V. CONCLUSIONS AND OUTLOOK 

We have investigated the ferroelectric trends of the 
incipient ferroelectric material Ti02 rutile by means of 
first-principles calculations. Similar to the well-known 
ferroelectric perovskites, rutile possesses anomalously 
large dynamic charges, which is one key factor for the 
stabilization of a ferroelectric phase. 

We could show that this strong dynamic charge trans- 
fer is mainly mediated by 7r-type 0-Ti hybrid orbit als. 
Furthermore, this dynamic charge transfer is consider- 
ably modified by changing the ratio of the equatorial and 
apical Ti-0 bond lengths. Nevertheless, a large dynamic 
charge transfer along the [001] direction and along the 
apical Ti-0 bonds is preserved under small variations of 
the lattice as they may appear at strained surfaces. In- 
deed, the Born charges in both directions are quite sta- 
ble at the (110) surface of rutile, even if one imposes an 
uniaxial tensile strain in the surface planes. Only small 
oscillations of the dipolar interaction appear with the sur- 
face distance. Here, the dipolar interaction along [001] is 
slightly enhanced within every second Ti-0 layer start- 
ing from the topmost layer, while the dipolar interaction 
along [110] is slightly larger in the odd layers. 

Although, the material possesses such a large dynamic 
charge transfer, no ferroelectric transition appears, as the 
short-range repulsion, which builds up during a polar 
displacement is too large for the undistorted material. 
In agreement with previous theoretical studies, we have 
shown that a ferroelectric transition of the bulk material 
can be enforced by lattice expansion or by uniaxial ten- 
sile strain. In both cases, the enlarged Ti-O distances 
reduce the short-range repulsion during a polar displace- 
ment. Additionally, we could show that a ferroelectric 
state polarized along [001] can be stabilized by all mod- 



ifications of the lattice which enlarge the nearest Ti-O 
distance, even if the apical Ti-0 distance shrinks during 
the lattice modification. For uniaxial strain along [001] 
and [1 10] stable ferroelectric states can be found and for 
both kinds of strain, the polar displacement along the 
strain direction is most favorable. 

Most notably, the stabilization of two ferroelectric 
states, one polarized along [001] and one along [110] can 
be transferred from the bulk material to the rutile (110) 
surface. We thus find large local dipole moments already 
for 2 % tensile strain in both [001] and [110] direction. 
Although the ferroelectric trends and the obtained local 
dipole moments may be slightly overestimated within our 
theoretical approach, we find a successive increase of the 
local dipole moments for increasing strain. Therefore, 
we predict that a ferroelectric state can be stabilized for 
experimentally accessible tensile strains, and that the re- 
sulting electric dipoles persist even in the atomic layers 
at or immediately below the surface. On the other hand, 
for compressive strain no stabilization of a ferroelectric 
phase is possible at the rutile surface, in agreement with 
the corresponding bulk material. 
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